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- - Abstract 

(N ■ 

The ANTARES telescope observes a full hemisphere of the sky all the time with a duty cycle close to 100%. This makes 
it well suited for an extensive observation of neutrinos produced in astrophysical transient sources. In the surrounding 
^ medium of blazars, i. e. active galactic nuclei with their jets pointing almost directly towards the observer, neutrinos may 
^ he produced together with gamma-rays by hadronic interactions, so a strong correlation between neutrinos and gamma- 
'rays emissions is expected. The time variability information of the studied source can be obtained by the gamma-ray 
^ light curves measured by the LAT instrument on-board the Fermi satellite. If the expected neutrino flux observation 
^^O is reduced to a narrow window around the assumed neutrino production period, the point-source sensitivity can be 

drastically improved. The ANTARES data collected in 2008 has been analysed looking for neutrinos detected in the 

high state period of ten bright and variable Fermi sources assuming that the neutrino emission follows the gamma-ray 
HH light curves. First results show a sensitivity improvement by a factor 2-3 with respect to a standard time-integrated 
. point source search. The analysis has been done with an unbinned method based on the minimization of a likelihood 
j-atio applied to data corresponding to a live time of 60 days. The width of the flaring periods ranges from 1 to 20 days. 
i" 1 Despite the fact that the most significant studied source is compatible with background fluctuations, recently detected 
O flares promise interesting future analyse. 
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1. Introduction 

Neutrino astronomy is an incipient field of observation 
of the universe. Since neutrinos do not interact electro- 
magnetically, but weakly only, they are not absorbed and 
point directly to their original sources, playing the role 
of unique messengers. However, due to their weak inter- 
action, neutrino telescopes are very low event rate experi- 
ments requiring high fluxes in order to detect a clear signal 
of a neutrino astrophysical source. Up to now no claim of 
an astronomy source of neutrinos further than the Sun has 
been made, with the exception of the well known SN1987A. 
Assuming a direct correlation between the neutrino emis- 
sion and the light emission of a source in an hadronic accel- 
eration scenario [l[ , the standard analysis can be improved 
by a factor 2-3 in the sensitivity if we limit this search to 
the periods of maximum probability of neutrino emission 
by the source. The time variability visible in the gamma- 
ray light curves measured by the Fermi-LAT Q can be 
used for this purpose. 

1.1. The ANTARES telescope 

The ANTARES detector is a neutrino telescope placed 
at the bottom of the Mediterranean Sea (42°48 N, 6°10 E), 
at a depth of 2475 m, connected by a submarine cable of 
42 km to the shore in Toulon (France). This cable connects 



through a junction box 12 lines which are separated by 60- 
70 m and vertically suspended by a buoy. Each line has 25 
floors spaced by 14.5 m, except for line 12 which has only 
20 floors. A floor consists of a triplet of optical modules 
(OMs) each one housing a photomultiplier (PMT) facing 
45° downwards. The full detector is a tri-dimensional ar- 
ray of 885 PMTs 3] 4] which was completed in 2008 when 
the last line was connected. 

Neutrinos are detected via the Cherenkov light induced 
by relativistic muons produced in the detector surround- 
ings by CC interactions of muon neutrinos with nuclei in 
water. The signals from the Cherenkov photons detected 
by the PMTs are digitized ('hits') [j| and sent to the shore 
station for reconstruction and physics analysis. 

2. Time-dependent analysis 

2.1. Analyzed data 

The analyzed data correspond to the period from Sep- 
tember 6th to December 31st, 2008 (54720-54831 modified 
Julian day). Some periods have been excluded from the 
selection in order to avoid conditions with high optical 
backgrounds due to the bioluminescence activity, leaving 
60.8 days of life-time. 

While most of the atmospheric muon background is 
suppressed by selecting only up-going events, a fraction 
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of mis-reconstructed events still remains which have to be 
rejected in a different way. In the track-reconstruction 
process an algorithm based on the maximization of a like- 
lihood is used. The likelihood function is built from the 
difference between the expected and the measured arrival 
times of the hits from the Cherenkov photons emitted 
along the muon track. Its maximization considers the 
Cherenkov photons that scatter in the water and the ad- 
ditional photons that are generated by secondary particles 
(e. g. electromagnetic showers created along the muon tra- 
jectory). A good measure of the track- fit quality is the 
value of the log- likelihood per degree of freedom (A), pro- 
viding a useful tool to reject the mis-reconstructed muons 
by applying a cut in that A value. This cut leaves the 
atmospheric neutrinos as the dominant background. An 
additional cut requiring the error estimated for the recon- 
structed muon track direction to be less than 1° is also 
applied in the selection of events. 

The angular resolution has been estimated by Monte 
Carlo simulations. The cumulative distribution of the an- 
gular difference between the reconstructed muon direction 
and the neutrino direction, with an assumed spectrum pro- 
portional to E y ~ 2 , where E v is the neutrino energy, is 
shown in Fig. [TJ For the studied period, the median reso- 
lution is estimated to be 0.4 ± 0.1 degrees. 

2.2. Time- dependent search algorithm 

An unbinned method based on a likelihood ratio max- 
imization has been used to perform this time-dependent 
point-source analysis. The data is parameterized as a two 
components mixture of signal and background. The goal 
is to determine, at a given point in the sky and at a given 
time, the relative contribution of each component and to 
calculate the probability to have a signal above a given 
background model. The likelihood ratio A is the ratio of 
the probability density for the hypothesis of background 
and signal (P s i g+ bkg) over the probability density of only 
background (if b kg) : 
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where n s i g and N are respectively the unknown number of 
signal events and the total number of events in the con- 
sidered data sample. P s - lg (ai,Si,ti) and Pbkg(cti, Si, U) are 
the probability density function (PDF) for signal and back- 
ground respectively, and Si is the declination of the studied 
source. For a given event i, ti and on represent the time 
of the event and the angular difference between the coor- 
dinate of this event and the studied source. 

P S i g (o!j, Si, ti) and Pykg(oii, Si, ti) are described by two 
components: 




Psi s (ai,Si,ti) = P s t(oH,Si) x P£™(U) 



(2) 



Figure 1: Cumulative distribution of the angle between the true 
Monte Carlo neutrino direction and the reconstructed muon direction 
for E~ 2 upgoing neutrino events selected for this analysis. 



Pbkg(ai,Si,ti) = P b kg (ati,5i) X Pbkg C (U) 

where P^ is the probability to have a signal event from 
the studied source, computed via Monte Carlo, P s '™ e is the 
probability of a neutrino event derived from the gamma- 
ray light curve correlation, and Pbk g ano - ^bkg C are the 
corresponding background PDFs. The shape of the time 
PDF for the signal event, P s '™ c , is extracted directly from 
the gamma-ray light curve assuming that it is proportional 
to the gamma-ray flux. For the signal event directional 
PDF, P s ii r , the one-dimensional point-spread function is 
used, which is the probability density of reconstructing 
an event at an angular distance on from the true source 
position. Then, the directional, P^ kg , and time PDF for 
the background, P b kg C , are derived from the data using re- 
spectively the observed declination distribution of selected 
events in the sample and the observed time distribution of 
all the reconstructed muons. The latter distribution (not 
normalized to 1) is shown in Fig. [51 where periods with- 
out data (due to detector maintenance, etc) or with very 
poor quality data (high bioluminescence activity or bad 
calibration) are shown as empty bins. 

The Adata value obtained from the data is then com- 
pared to the distribution of A given by the null hypothesis 
(n s ; g = 0). The comparison of Adata with the background 
only A distribution is used to reject the null hypothesis, 
where the confidence level corresponds to the fraction of 
the scrambled trials above Adata (p- value). The discovery 
potential is then defined as the average number of signal 
events required to achieve a p- value lower than the equiv- 
alent to 5<T in 50% of trials. 

The average number of events required for a 5cr dis- 
covery (50% C.L.) produced in one source located at a 
declination of -40° as a function of the total width of the 
flare period is shown in Fig. [31 A comparision with the 
numbers obtained with a standard time-integrated point 
source search, shows an improvement of the discovery po- 
tential by about a factor 2-3. 
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Figure 2: Time distribution of the reconstructed events. Black: dis- 
tribution for all reconstructed events. Red filled: distribution of 
selected upgoing events (A > —5.4 and < 1°). If 0, there are no 
data available (i. e. detector in maintenance) or the data have a very 
poor quality (high bioluminescence or bad calibration). 
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Figure 4: Gamma-ray light curve (red dots) of the blazar 3C454.3 
measured by the LAT instrument on board the Fermi satellite above 
100 MeV for almost two years of data. Blue histogram: high state 
periods. Green line and dots: baseline and significant dots above 
this baseline used for the determination of the flare periods. 
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Figure 3: Average number of events required for a 5cr discovery (50% 
C.L.) of a source located at a declination of -40° as a function of the 
width of the flare period (solid line), a simple Heavyside function. 
This is compared to the number required for a time-integrated search 
(dashed line). The average typical length of the studied source flares 
goes from 1 up to 100 days. 

2.3. Search for neutrino emission from gamma-ray flares 

The sources for which this time-dependent analysis has 
been applied have been selected from the Fermi blazar 
sources reported in the first year Fermi LAT cataloguelg] 
and in the LBAS catalogue (LAT Bright AGN sample 0). 
The selection has been done based on their variability, 
brightness and visibility: sources visible to ANTARES 
with a significant time-variability in the studied period 
and a flux in the high state greater than 20 x 10~ 8 pho- 
tons cm" 2 s _1 above 300 MeV in the averaged 1 day-binned 
light curve were selected. The final selection list of sources 
includes four BLLacs and six flat spectrum radio quasars 
(FSRQ), listed in Table [TJ Table Q] lists their fluxes, visi- 
bility and live times. 

The light curves used for P s '™° in j2} are those pub- 
lished on the Fermi web page for monitored sources Q 
(one-day-binned time evolution of the average gamma-ray 
flux above a threshold of 100 MeV). Figure 0] shows the 
high state periods (blue histogram) obtained for the Fermi 



LAT gamma-ray light curve of 3C454 for almost two years 
of data. These states arc defined using a simple and robust 
method based on three steps: 

• First, the baseline and its error is determined with 
an iterative linear and gaussian fits. After each fit, 
the points where the flux is above the fitted base line 
plus one sigma are suppressed. This is done 3 times. 

• Then all points (green dots in Fig. 0]) where the mea- 
sured flux minus its error is above the baseline plus 
two times the baseline sigma and at the same time 
the measured flux value is above the baseline plus 
three times the baseline sigma, are used as priors 
from which the flares are defined. 

• Then, for each selected prior, the adjacent points 
for which the emission is compatible with a flare are 
added, i. e. the adjacent points which have their flux 
minus their error above the baseline plus one sigma. 
Finally, an additional delay of 0.5 days is added be- 
fore and after the flare in order to take into account 
the 1-day binning of the light curve. Hence, a flare 
has a minimum width of 2 days. 

Assuming that the neutrino emission follows that in gamma- 
rays, the signal time PDF is simply this de-noised light 
curve after normalization. 

3. Results 

The most significant source is 3C279, which has a pre- 
trial p-value of 1.03%. The unbinned method finds one 
high-energy neutrino event located 0.56° from the source 
location during a large flare in November 2008 (Fig. @|. 
This event has been reconstructed with 89 hits spread on 
10 lines with a track fit quality A = —4.4 and an error es- 
timate of j3 = 0.3°. The post-trial probability is computed 
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Source 


OFGL name 


Class 


Redshift 


-F300 


Visibility 


Live Time 


N(5a) 


A obs 


Fluence 


PKS 0208-512 


J0210.8-5100 


FSRQ 


1.003 


4.43 


1.00 


8.8 


4.5 





2.8 


AO 0235+164 


J0238.6+1636 


BLLac 


0.940 


13.19 


0.51 


24.5 


4.3 





18.7 


PKS 0454-234 


J457.1-2325 


FSRQ 


1.003 


13.56 


0.63 


6.0 


3.3 





2.9 


OJ 287 


J0855.4+2009 


BLLac 


0.306 


2.48 


0.39 


4.3 


3.9 





3.4 


WComae 


J1221.7+28.14 


BLLac 


0.102 


2.58 


0.33 


3.9 


3.8 





3.6 


3C 273 


J1229. 1+0202 


FSRQ 


0.158 


8.68 


0.49 


2.4 


2.5 





1.1 


3C 279 


J1256.1-0548 


FSRQ 


0.536 


15.69 


0.53 


13.8 


5.0 


1 


2.8 


PKS 1510-089 


J1512. 7-0905 


FSRQ 


0.360 


28.67 


0.55 


4.9 


3.8 





2.8 


3C 454.3 


J2254.0+1609 


FSRQ 


0.859 


24.58 


0.41 


30.8 


4.4 





23.5 


PKS 2155-304 


J2158. 8-3014 


BLLac 


0.116 


7.89 


0.68 


3.1 


3.7 





1.6 



Tabic 1: List of the bright, variable Fermi blazars selected for this analysis. -F300 is the gamma-ray flux above 300 MeV in 10 — 8 photons cm 2 s _1 . 
Live Time is the effective time of observation of the source in days. N(5cr) are the number of neutrino events needed to be detected in 
ANTARES with 5cr, while N is the number of events observed in the source's direction and flare time window. Fluence is the upper limit 
(90% C.L.) on the neutrino fluence in GeV cm -2 , calculated according to the classical (frequentist) method for upper limits 0. 
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Figure 5: Gamma-ray light curve (black dots) of the blazar 3C279 
measured by the LAT instrument on board the Fermi satellite above 
100 MeV. Blue histogram: high state periods. Green dashed line: fit 
of a baseline. Red histogram: time of the ANTARES neutrino event 
in coincidence with 3C279. 

taking into account the ten searches. The final probabil- 
ity of 10% is compatible with a background fluctuation. 
Other source results are sumarized in Table [TJ 

4. Conclusion 

The first time-dependent search for cosmic neutrinos 
in ANTARES has been presented. It has used data taken 
with the full 12-line ANTARES detector during the last 
four months of 2008. Time-dependent searches are signifi- 
cantly more sensitive than a standard point-source search 
for a variable source. This search has been applied to ten 
bright and variable Fermi LAT blazars. The most signifi- 
cant observation of a flare is 3C279 with a p- value of about 
10% after trials for which one neutrino event has been 
detected in time and space coincident with the gamma- 
ray emission. Limits on the neutrino fluence have been 
obtained for the ten selected sources. The most recent 
measurements of Fermi in 2009-11 show very large flares 



yielding a more promising search for neutrinos 10], and 
ongoing analyses are being performed together with im- 
provements in the denoising of the light curves. 
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